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Abstract 

We address in this paper the approximation problem of distributed delays. 
Such elements are convolution operators with kernel having bounded support, 
and appear in the control of time-delay systems. From the rich literature on 
this topic, we propose a general methodology to achieve such an approxima- 
tion. For this, we enclose the approximation problem in the graph topology, 
and work with the norm defined over the convolution Banach algebra. The 
class of rational approximates is described, and a constructive approxima- 
tion is proposed. Analysis in time and frequency domains is provided. This 
methodology is illustrated on the stabilization control problem, for which 
simulations results show the effectiveness of the proposed methodology. 
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1. Introduction 

The interest for the use of distributed delays in the stabilization of time- 
delay systems appears in the pioneering work of Olbrot [lj]. To generalize 
algebraic methods issued from linear systems in finite dimensional spaces to 



* Corresponding author 
Email addresses: hao.lu@insa-lyon.fr (Hao Lu), 
michael.di-loreto@insa-lyon.fr (Michael Di Loreto), 

damien.eberard@insa-lyon.fr (Damien Eberard), jean-pierre.simon@insa-lyon.fr 

(Jean-Pierre Simon) 

1 The authors acknowledge the financial support of the CSC Scholarship from Xi'an 
Northwestern Polytechnical University (China), and the scholarship program with INSA- 
Lyon (France). 



1 



time-delay systems, Kamen et al. [2] first introduce a general mathematical 
setting for the control, and in particular for the stabilization, of time-delay 
systems. This mathematical framework was formalized in [3] by the intro- 
duction of the Bezout ring of pseudopolynomials, or independently by j3] 
with a behavioral approach. In these works, distributed delays appear in 
the stabilization problem. More generally, they are at the core in spectrum 
assignment by feedback j^] . Finite spectrum assignment generalizes the prin- 
ciple of Smith's predictor for dead-time systems jsj to more general systems 
with delays, which can be stable or unstable. Distributed delays appear also 
in the characterization of equivalence transformations 0. Robustness issue 
and optimization jsj, robustness for input-delay systems jsj or finite time 
control lOj are other topics where distributed delays play a central role. 



A distributed delay is a linear input-output convolution operator of the form 



y(t) = (f*u)(t) 



f{r)u(t-T)dr 



(1) 



where i? is a strictly positive and finite real, and kernel /(•) is a continuous 
function with support [0,1?]. Numerical implementation of distributed delay 
was early investigated, to propose effective algorithms for control. Such an 
implementation was the starting point of a widely research activity. A first 
proposition for approximation with finite dimensional systems was proposed 
in [11] . Reduction and approximation of delay systems, involving lumped 



delays, were also investigated in [12]. In the work of [13] . the authors pro- 
pose a numerical integral approximation to realize an operator like in ([[]). 
Such an approximation writes as a sum of lumped delayed distributions, 
and unfortunately introduces additional closed-loop poles, and also insta- 
bility phenomena 



See, e.g. [14], [15 



and references therein. To overcome 
it was outlined that 
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this problem, various solutions were proposed. In 
such an approximation has a poor accuracy and an high sensitivity for high 
frequencies. Hence, the author proposed to add a low-pass filter in the inte- 



gral approximation. Such a solution was also proposed indepedently by [17 
Further implementation improvements were proposed in 18], with rational 
approximation and extension of bilinear transformations. These last papers 
give positive answers to the open problem of general approximation of dis- 



tributed delays, outlined in [19 



We propose in this paper a general methodology for numerical implemen- 
tation of distributed delays. The numerical implementation of an operator 
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exhibit two sides. The first one involves time discretization of the input- 
output behavior of the operator. Any distributed delay is a BIBO-stable 
operator. It is easy to verify that an equivalent discrete time system can be 
obtained by usual sampling tools, and that this system can be put into a sum 
of causal lumped distributions. Taking an appropriate sampling period, this 
equivalent discrete time system can always be made BIBO-stable. Hence, this 



first part presents no difficulty We refer for instance to [20| and in references 
therein for more details on this part. The second part involves the approxi- 
mation problem in continuous time of such an operator. This is this part we 
address in this paper. A continuous time approximation need to reproduce 
with high fidelity the internal dynamics of this operator, for large classes of 
input signals, but also to generate an arbitrarily close input-output behav- 
ior to the original one. For linear systems governed by a convolution, these 
notions are equivalent to impose simultaneous time and frequency approxi- 
mations. For this, we explicitly separate two notions, namely input-output 
approximation and kernel approximation. We will explain why input-output 
approximation is not suitable for approximation of distributed delays, and we 
will then focus ourselves on kernel approximation. Kernel approximation can 
be realized in many different ways. Among most used classes of operators for 
approximation, we can cite polynomials, rational fractions, or exponentials. 



See, e.g. [2jj, [22|, or [23 



With the objective to substitute the distributed delay by a more tractable 
system, highlighting rational assumptions, we propose two classes to realize 
approximation, namely lumped systems and a subclass of distributed de- 
lays. With previous objectives, we enclose the approximation problem into 
the Wiener algebra of BIBO-stable systems, using the graph topology. This 
corresponds to the weakest topology where feedback is a robust property. 
Moreover, for stable systems, graph topology and norm topology being the 
same, we work on norm convergence over this algebra, which is a Banach 



algebra. This general framework was used first in [24[ for approximation of 
distributed parameter systems by lumped systems. Roughly speaking, work- 
ing over this algebra, an approximation comes down to realize approximation 



of the kernel over the Banach algebra Jzfi(R + ). This idea also grew in [25 
for the approximation of lumped delayed distributions which appear in op- 
timal control. Here, we propose an extension of the classes of approximates, 
and we show that working in this general setting yields to approximate a 
distributed delay in both time and frequency domains, for large classes of 
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input signals. We also propose to highlight the propositions made in [16 
or 17| . and to bring a mathematical foundation for this solution. 



The paper is organized as follows. In Section [2j we define and fully char- 
acterize the main properties, in both time and frequency domains, of dis- 
tributed delays. We characterize in particular a general decomposition of 
distributed delays on the so-called elementary distributed delays. In Sec- 
tion El we explicit our approximation problem, and solve it. This section 
starts with general comments on approximation, where graph topology is 
briefly recalled, and where we show that any distributed delay can be ap- 
proximated over this topology by lumped systems. Then, we move to another 
class of approximates, using a subclass of distributed delays, easy to imple- 
ment with stability. We show a density property of this subclass, and analyze 
the convergence of the approximation which is proposed. Section H] relates 
the properties of such an approximation. We outline a proposal for a con- 
structive approximation, and analyze time and frequency properties. A few 
simulation show the effectiveness of the method, on the stabilization control 
problem. 

2. Convolution operators and distributed delays 

2.1. Convolution algebra 

An input-output causal convolution system is a dynamical system de- 
scribed by an equation of the form 



where ?/(•), u(-) and /(■) are said to be the output, input and kernel of the 
map, respectively. Convolution systems are naturally defined over a commu- 
tative algebra since they are closed under addition, multiplication and scalar 
multiplication, operations that correspond to arbitrary series and parallel in- 
terconnections of such systems. For a normed algebra, it is quite interesting 
to obtain some closure properties for convergence. Hence, we prefer to work 
on a Banach algebra, that is an algebra for which any convergent sequence 
of elements in the algebra has a limit in the algebra, and for which the norm 
has the multiplicativity property. It will be assumed that all linear spaces 
and algebras are over the complex field. A general algebra of distributions 
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including a wide class of convolution systems is given by the so-called Callier 
and Desoer algebra, denoted stf 26[]. We say that / G if 



/<») = { £ (t)+/ - (<) " HI (3) 

where the complex-valued function /„(■) G Jzfi(R + ), that is / a is a complex 
valued function, locally integrable on R, + , and such that J °° \f a (t)\dt < oo. 
The complex-valued distribution f pa stands for the purely atomic part and 
writes 

oo 

f P a(t) = Y t US(t-t n ), (4) 
n=0 

with f n E C, n = 0,1, = t < t\ < t 2 <■■■, 5(t — t n ) denotes the Dirac 
delta distribution centered in t n , and J2 n >o \fn\ < oo. As shown in Desoer 



and Vidyasagar |27j, it is well known that srf is a commutative convolution 



Banach algebra with norm defined by 

oo 



si 

n=0 



and with unit element the Dirac delta distribution 5. Denoting / the Laplace 
transform of /, srf denotes the set of Laplace transforms of elements in srf . 
The set srf is also a commutative Banach algebra with unit element under 
pointwise addition and multiplication, for the norm 



The algebra srf gives a general mathematical framework for the analysis of 
distributed delays, as it will be explained in the next subsection. Before this, 



let us recall the concept of Bounded Input-Bounded Output stability [27 



Definition 1. A convolution system in the form (T5]) is said to be BIBO 
stable if / G srf . 

More generally, we may be interested by Jz^,-stability, for 1 < p < oo. We say 
that (j2J) is Jz^,-stable if for any u(-) G =^(R + ), that is u(-) locally integrable 
over R + and 

i/p 



\uy p = Ht)\ p d?j < oo, 
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the output y(-) is also in Jz^,(R + ). Since 

IMk = 11/ * «IU P < ll/IUIkll^,, (6) 



^-stability is equivalent to BIBO stability [27 



2.2. Distributed delays 

Let \ a ,b = [a,b] be the bounded closed interval in R + , for some reals a 
and b, < a < b. Notations I 0)OO or R + stand for [0, oof. We define J(f{\ a ^) 
as the set of complex valued functions g(-) of the form 

g( t) = {^M tel* (7) 

w [0, elsewhere v ' 

where 

i>0 j>0 

for some cy and Aj in C, and the sums are finite. In other words, gi a b is 
a finite linear combination of exponential-polynomials type functions, and 
it is in particular a continuous function. For any real valued function in 
Jf(I at b), if some A; 6 C appears in the sum, then so does its conjugate Aj, 
and the associated coefficients Cjj are complex conjugates. Hence, any real 
valued function in J{f{l a ^) is a function generated by real linear combinations 
of Pe ait , Pe ait sxn(f3kt) and Pe ait cos((3kt), for some real numbers <Tj, 0k, the 
sums being finite. The formal definition of distributed delay is made below. 

Definition 2. A distributed delay is a causal convolution system, with ker- 
nel / in ^ 2 ), for some bounded real numbers < $1 < 

In other words, a distributed delay can be written like an input-output con- 
volution operator of the form 

y(t) = (f * u)(t) = h^ 2 (r) U (t-r)dT, (9) 

with notations introduced in ([8]). The set of distributed delays, denoted by 
^, is a ring. Obviously, for real valued signals, the kernel /(■) will be a 
real valued function in Jtr(l# lt # 2 ). In the previous definition, we restrict our- 
selves to define a distributed delay like a convolution operator with kernel in 
^(Itfii? 2 )- This restriction is not conservative as we will show in Section 1X21 



6 



Actually, all distributed delays which appeared in the literature are particu- 
lar cases of this definition. This definition is based on a rational construction, 
as this appears explicitly using Laplace transforms. Any distributed delay £f 
admits a Laplace transform, corresponding to the finite Laplace transform of 
its kernel / £ JtT(I$ u .& 2 ), 

y(s) = f(s)u(s), f(s) = I' (r) e — dr, (10) 

where / £ £f is an entire function, i.e. holomorphic on the whole complex 
plane. The notion of elementary distributed delay will greatly simplify the 
approximation problem. Let us define the complex valued function 6\(-) £ 
J^(I 0) #), for some A £ C and $ > 0, by 



e ^>-S 0, elsewhere (U) 



and its Laplace transform 



0x(s) = e —r— , (12) 

s — A 



which is an entire function even in s = A where 9\(\) = In other words, A is 
a removable singularity, and consequently 9\{s) has no pole. The distributed 
delay whose kernel is 6\ is called an elementary distributed delay. The kih. 
derivative * \s) of 9\{s) yields 

9[ k \s)= [\-r) k e-^dr, (13) 
Jo 

which is still in £f , and corresponds to the Laplace transform of the function 
9\{t) = (—t) k e xt for t £ [0,19], and elsewhere. From previous definitions, 



we can state the following lemma, which also appeared in [28( using 2D- 
polynomials, and which will play a central role for approximation. 

Lemma 1. Any element in^ can be decomposed into a finite sum of Laplace 
transform of elementary distributed delays and its successive derivatives. 

PROOF. Take any element in . Its kernel g(-) lies in Jf{J.^ 1 ^ 2 ), and writes 
as in ([8]). By time translation corresponding to the lumped delay it is 
readily a linear finite combination of elementary distributed delays 9\(t) and 
of the functions 9\{t), as defined in ffTTl) and ([TBI , with •§ = $ 2 — 1=1 
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In other terms, previous result tells us that for any g G there exist complex 
polynomials g^ G C[e _1?s ] with respect to the variable e~® s and Aj G C, in 
finite number, such that 

^) = E^( e ~" s )C } (^ (14) 

i,k 

where successive derivatives are iteratively computed by 

0£\s) = (-1)**!— b^€r, (15) 



and 0^(s) are still entire functions, since 0^(Ai) = ^-^n + , for any k > 0. 
From (I14p and taking into account that practical distributed delay is a real 
valued operator, we see that, grouping the terms via least common multiple, 
any element g G can be put into a fraction 

9{s) = -d&- ' (16) 

where n(s,e _,3s ) G R[s,e _1?s ] is a real quasipolynomial with respect to the 
algebraically independent variables s and e~^ s , and d(s) G R[s]. Any element 
in the right hand side of (|14p being an entire function, g(s) is also an entire 
function. Hence, any zero of d(s) is also a zero of n(s, e~® s ). Furthermore, the 
degree with respect to s of 0^\ s ) i n (EG]) and f lT5|) . for any h, being strictly 
negative, the degree with respect to s of g(s) satisfies deg s n < deg s d, so 
that its Laplace transform is strictly proper. As in [3|, the ring ^ is the 
ring of those Laplace transforms of distributed delays that are rational in 
the variable s and e _1?s , which are entire and strictly proper (with respect 
to s). This result comes from the assumption on rational kernels in (jS]). 
The kernel of any element in is obviously in j£fi(R + ), so that from (J6]) 
and the definition of =sz/, any distributed delay is BIBO stable. Convolution 
of two kernels with finite support yields another kernel with finite support, 
and since exponential-polynomials type functions are closed under classical 
product, we have in fact that £f is a normed subalgebra of srf for the || ■ ||^- 
norm. From these definitions and properties, we are now able to formulate 
the approximation problem of distributed delays. 
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3. Approximation of distributed delays 

For this purpose, we define two subspaces in direct sum in J^(I^ lt ^), 
denoted respectively JtgQL^^) and J^I^,^), consisting of linear combi- 
nations of exponential-polynomials type functions on some finite interval as 
in OH]), with Re Aj < and Re Aj > 0, respectively, for all i > 0. 
From Lemma [Hand CHI) , we know that any element g 6 J^(I# lj ^ 2 ) is a linear 
combination of elements 9\. in J(f s {\^ 1 ^ 2 ) 1 for which ReAj < 0. Hence, any 
element in J^(I^ lj ^ 2 ) can be numerically implemented with stability using el- 
ementary blocks, as illustrated in Fig. [TJfor 9\(t). Note that in practice, since 
A G C, we should implement it using a real decomposition. Approximation 
for the subclass J^il^^) of distributed delays is not required, since this 
implementation is realized with stability, and requires, from an implementa- 
tion point of view, only two pointwise delays, namely and However, this 
is no more true for elements in J(f u (I^ 1 ^ 2 ), where numerical realization for 
9\{t) as in Fig. [T] yields to an unstable system. Note that another realization 
for 9\{s) may be 

x(t) = Xx(t) + u(t) - e M u(t - 0), 

which is however still numerically unstable for Re A > 0. Intuitively, this 
instability is a consequence of a non exact numerical cancellation of the 
removable singularity s = A of 9\(s). Therefore, we will focus our attention 
on distributed delays whose kernels are in J^ U (I &1 ^ 2 ). 
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Figure 1: Realization of an element 9\(s) G J^(Iq,#)) with Re A < 0. 

The question can now be stated as follows: How to define a continuous time 
approximation of a distributed delay whose kernel lies in J^ U {1^ 1 ^ 2 )1 We 
give complete answers to this question in the next two subsections. We start 
by some general remarks on approximation over s$ ', and hence give answers 
to the previous question. 
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3.1. General comments on approximation 

Any distributed delay with kernel / G ^(I^ l5 ^ 2 ) is a convolution opera- 
tor. Its J^L-induced norm is 



si 



From ([6]), we see that this norm is actually an upper bound of all its induced 
££ v - norms, for 1 < p < oo. A natural metric is then obtained from the 
■cZ-norm, which in turn can be enclosed in the graph topology. Consider a 
distribution p with a coprime factorization (n, d) in A, that is p = hd" 1 , and 
such that there exist x and y in srf satisfying 

n*x + d*y = 5, (17) 

or in the Laplace domain 

nx + dy = 1, Vs G C, Res > 0. (18) 

A neighbourhood of p in the graph topology is the set of all plants of the 
form n^d^ 1 where (n^, d&) belongs to some ball in srf centered at (n, d). If p 
is itself in for any positive number r, a neighbourhood of p in the graph 
topology writes 

&(p,r) = {p A E £/ : \\p- Pa\\^ <r}. 



In [29|, Ch. 7], it was shown that the graph topology is the weakest topol- 
ogy in which feedback stability is a robust property. Graph topology being 
metrizable and hence first- count able, for BIBO-stable systems, norm topol- 
ogy and graph topology are the same. Let e be a given suitable small positive 
number. Hence, we say that a distribution p £ in A is an approximation of p 
in the graph topology if the element p e is close to p, that is p £ G £$(p, s), or 
equivalently \\p — p e \\&? < e. This property can be related to convergence in 
but we need first to define the class of operators that will approximate p 
in the graph topology. The most commonly used class is the so-called class 
of lumped systems, which consists in convolution operators whose kernels lie 
in the set £}, defined by 

■2 = {g G : g{t) = g 5{t) + g a (t), g G C, g a G J^(I 0) oc)} • (19) 



10 



From 24J, we know that any plant p with coprime factorization (n,d) over 
s$ can be approximated by a lumped plant in the graph topology if and only 
if there exists a real constant matrix M of rank 1 such that 



MS pa = 0, where S pa 



dpa 



(20) 



and {')pa denotes the purely atomic part. Distributed delays having strictly 
proper transfer functions, we are interested for approximation with the class 
of strictly proper lumped systems, denoted £ s , set of elements in J^(I 0jO o)- 
Such elements satisfy ordinary differential equations, and are easily simu- 
lated. From j23[, we know that the closure of £} s is j£fi(R + ), since every 
function in = Sf 1 (R + ) can be approximated by a sum of exponentials. This 
yields the following result, that can be seen as a particular case of the above 
given condition (|20|) . 

Theorem 1. Any plant in =Sfi(R + ) can be approximated by a lumped system 
in the graph topology. 

Proof. Let / be an element in «£?i(R + ). A coprime factorization of / over 
s$ is n = / and d = 5. Since n pa = 0, any real constant matrix M = [0 *], 



(*) denoting an arbitrary real number, satisfies (|20|) . Hence, from 24j, Th. 



4.1], / can be approximated by a lumped system in the graph topology. □ 



Note that if the element / in Jzfi(R + ) is strictly proper, then an approxi- 
mation in the graph topology realized by a lumped system will be in fact 
in the set £} s . Any convolution system whose kernel lies in j2fi(R + ) can be 
approximated by a lumped system. Roughly speaking, since any distributed 
delay is a strictly proper fraction and is BIBO-stable, we know that it can 
be approximated in the graph topology by a strictly proper lumped system. 



Corollary 1. Any distributed delay can be approximated by a lumped system 
in the graph topology. 

Others approximation classes can be used to numerically realize a distributed 
delay. In the literature, appears a method based on a numerical integral ap- 
proximation, but such a method was shown to be at the core of various unde- 
sired results, like instability or numerical sensitivity. Let us briefly make some 
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considerations on such approximations, introducing what we called input- 
output approximation, and let us show why such approximations are not 
suitable for our problem. For any distributed delay y(t) — (f * u)(t), we say 
that y app is an ££ p input-output approximation of y, if for some arbitrarily 
e > 0, \\y — 2/app 1 1 < £ holds. For instance, consider the Jzfoo input-output 
approximation based on integral approximation, like Newton-Cotes quadra- 
ture methods, which leads to an approximation y app (t) of the form 

g 

ya PP (*) = ^2 c kf( T k)u(t ~ r -0' ( 21 ) 

where reals Ck and depend on the applied method [3(3, Ch. 4]. Various 
problems arise with this kind of approximation. First note that such an ap- 
proximation is made for a given u(-). Hence, if such an input changes, the 
properties of this approximation failed, in general. Note also that, as men- 



tioned by 16|, the Laplace transform of (j2ip yields to an input-output ap- 
proximation transfer function that is proper, since it contains only pointwise 
delays. From (j!4p and ffl6|) . a distributed having a strictly proper transfer 
function, this approximation has a poor accuracy for high frequencies, which 
will have as consequence a high sensitivity to high-frequency plant uncertain- 
ties, in particular for the closed-loop system. This fact is illustrated in Fig. [21 
where a frequency diagram of 0\(s) and an input-output approximation is 
plotted. This negative result is directly interpretable in the graph topology 
by the following result. 

Theorem 2. Any distributed delay can not be approximated in the graph 
topology by a purely atomic distribution. 

Proof. Consider an element in with kernel g. Then g e Jzfi(R + ). Con- 
sider an arbitrary purely atomic distribution with pointwise delays of the 
form 

i>0 

where J2 i>0 \pi\ < oo. Such a distribution corresponds to the kernel of the 
approximation using numerical integral approximation like in (j21j) . Then 



\g-p\U = IMU + \ pi 



i>0 
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Figure 2: Bode gain diagram of 9i(s) and an approximation based on integral approxima- 
tion (trapezoidal Newton-Cotes method with 30 nodes). 

Hence \\g — p\\j# > ||g||j?i, and consequently p can not approximate g in the 
graph topology. □ 

It is clear that input-output approximation is not suitable for numerical 



implementation. This was illustrated in [13[, [14J, or [16j. Hence we focus 
our attention to the notion of approximation defined in the graph topology. 
This approximation is equivalent to realize an approximation of the kernel of 
the convolution, so we may also call it kernel approximation, to discern with 
respect to input-output approximation. 

3.2. Kernel approximations 

Approximation in the graph topology for an element / in Jtf s {\$ 1 ,Q 2 ) re- 
quires to find, for some given e > 0, an element / app , defined in general in 
such that / app G &(f,e). Element / being in particular in J>fi(]R + ), 
adding purely atomic part in approximation can not improve approximation, 
so we realize it over the Banach algebra Jzfi(R + ). Indeed, if an approximation 
writes 

f .._ fapp i xapp 
J app J a ~ J pa J 

where f% pp and fp PP stand for the atomic and purely atomic parts of / app , 
respectively, the following decomposition holds 

II/ — /applliz/ = ||/ — fa PP \\^ + ll/pa P H^- 
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The purely atomic part of the approximation is in an independent sum, and 
can then be reduced to zero to reduce the approximation error. We have 
seen that an approximation can be obtained by lumped systems, with ratio- 
nal and stable transfer functions. Since we want to obtain an approximation 
which is global with respect to the time, we are interested in uniform con- 
vergent sequence of functions in «if 1 (R + ) to /. Consider a sequence /„ of 
functions which uniformly converges to /, that is, for any e > 0, there exists 
n in IN such that for any k > n, \\f — f K \\^ < e. So, for any k > n, f K is a 
kernel approximation of / in the graph topology. We propose next a method 
to approximate an element in J£u(I# 1} & 2 ) by elements in J^I^^J. Since 
these last elements are easily implementable, this will give us an effective 
approximation. In other words, we approximate a distributed delay by an- 
other distributed delays. We start by considering the case of an elementary 
distributed delay. 

Lemma 2. Any distributed delay with kernel 0\(-) in J^(I 0)1 ?) can be approx- 
imated by distributed delays with kernels in J^(Io^) for the graph topology. 

Proof. Let 9\(-) be the kernel of an elementary distributed delay, with 
Re A > 0. Let \l £ C such that Re \i = a < 0. We define the transform 

A (a) = (as)- 1 ^(-a" 1 In s), e~* a < s < 1. (22) 

The function 0\(-) is continuous over [e _1?a , 1], and we have 

||0a||jz?i = II^aII^- (23) 

By Miintz-Szasz theorem, we see that Q\(-) can be approximated with respect 
to „2i-norm, as closed as desired, by a function of the form 

e A ,app(a) = (as)- 1 *A,a P p(s), (24) 

where ^x tapp is a polynomial in s. Since norms are preserved in f !2"3"j) . this 
comes down to approximate 6\ as closely as we please, using inverse transform 
of (|22|) . by a sum of exponentials of the form 

^A,app(t) = ^A,app(e~ at )> 

which is clearly an element in Jf s {J.Q^). □ 
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In other words, an approximation of 9\ in stf can be a function of the form 

n 

^A,n(*) = X^V»M*)' ( 25 ) 
i=l 

with ctj G C, Recti < some arbitrarily complex numbers, and 7^ some 
suitable constants. The order of the approximation n describes the number 
of parallel distributed elements to be added, to get the desired approxima- 
tion accuracy. Such an approximation does not increase the degree of the 
elementary distributed delay, since it is a sum of distributed delays which 
are strictly proper and have same degree than 9\(s). From the synthesis 
of an approximation for the elementary distributed delay, we show in the 
next two results that it allows to construct an explicit approximation for any 
distributed delay. 

Lemma 3. Let ipx,n be an element in J^(Io,#) H &(9\, e), for a given e > 0, 
and k in IN. Then G J^(I ,tf) n SS{d\,4 k e), where 

^ {t) = { i-t) k 1>X,n(t) , te[O,0] 

A,n \ 5 elsewhere 
Proof. Let e > 0, and take ipx,n an element in J^(I 0) #) fl &(0\,e), that is 

n 

l/>\,n(t) = ^li, n 9 ai {t), 
i=l 

with 7 ij7l in C, and \\6\ — ip\ tn \\,^ < e. By Laplace transform and fcth order 
differentiation, we get 

$S(«) = £7«,» *<!?(«)■ 
i=i 

In the time domain, this last identity corresponds to 

n 

<n(o = £7iX(o = (-o fc iM*)- 



i=l 



Clearly G J^ S (I ^). Hence, it remains to show that ip^ n G g§{6\,f> k e). 
But ip\ jn G 3S(Q\,e), so we have 



i9 



l|0A-<nlU = / |H)*(flA(0-^A,n(*))|d* 
< ^||^-^A,n|U<^e, 
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which completes the proof. 

For any distributed delay, we have the following general result. 
Theorem 3. Let g e £f be an arbitrary distributed delay, of the form 

= £^e-^%), 



□ 



9{s) 



i.k 



and let £j be given positive real numbers. For any elements ip\ i)n in J^ S (I ^) fl 
.^{Ox^Si) , we define 



2app(s) = ^Me-^2; n (s) 



Then g app lies in J^(Iq^) D &(g, e), where e = M ■ maxj £j /or some positive 
constant M. 



Proof. Let g be the kernel of a distributed delay. According to (fl4j) . we 
decompose it as a linear combination of elementary distributed delays and 
their successive derivatives. From the definition of Qa P p 5 which is clearly in 
, we have 



\9 ~~ 9a,pp\\if 



E^( e "*')(^( s )-%w) 



t.fc 



Since, for any z, are m £${9\v e i)i using Lemma EJ we get 



Denote the positive bounded constant M = J^ifc^ll&fcC 61 Hence 

\\g - OapplU <e = M ■ max e h 

i 

that is g app lies in £8(g,e). 



□ 
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By suitable choices for e i} the upper bound e can be reduced arbitrarily, so 
that we can find g app as close as we please of g. Note that we particularize 
the proofs of previous results with approximations with kernels in J^(Io^), 
but they can be trivially extended to others approximations over stf ' . Previ- 
ous results state that from the approximation of the elementary distributed 
delay, we can realize a kernel approximation in the graph topology for any 
distributed delay in Previous results state that, for any / G ^(I^^) 
and e > 0, there exists an element / app G ^(I^ l5 ^ 2 ) such that 



Said differently, we have the following corollary. 

Corollary 2. The set Jfc("l& u # 3 ) is dense in Jf(J.0 lt # 2 ) for the graph topol- 
ogy- 



In the proof of Lemma El convergence of polynomial approximation is uni- 
form, so we get here a uniform convergence of this approximation. The 
assumption that kernel of distributed delay lie in J^(I^ U ^ 2 ) is not restrictive. 
Indeed, if this is not the case, any continuous function g in [^i,^] can be 
approximated with respect to Jzfx-norm as close as we please by a function in 
=-^(Itfi,i? 2 ) @1 ■ Previous results addressed the approximation of distributed 
delays using kernels in J^(l^ u ^ 2 ), which is based on uniform convergence of 
polynomials to any continuous function over I# u a . Other approximation can 
be proposed. Indeed, any element / in Jif(TL0 u 2 ) can also be uniformly ap- 
proximated by stepwise continuous functions in J£?i(Io lCO ). This can be seen 
as a particular case of the previous approximation. For this, consider 



app 



(26) 



Proof. Obvious from Lemma [T]and Theorem [3j 



□ 



n 




(27) 



i=0 



where ^ = — , and ip(-) is a function in J^(I 0jO o)- Then 






n-l 



9x(t)-J2^(t-U)\dt+ / |0A,app(t)|dfc 
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From Lemma [21 we see that approximation (j27p uniformly converges to 6\(-) 
in the graph topology by a suitable choice of coefficients 7^. Taking ip(t) = 
e~ at h(t), where a > and h(-) stands for the Heaviside function, we obtain 



as a particular case the result obtained in [16| , and separately in [17| , where 
a low pass filter is added in the integral approximation with lumped delayed 
distributions. Indeed, Laplace transform of ( 12"T|) is of the required form 



1 " 



s + a 

for a > 0. Coefficients 7^ can be obtained for instance from some numerical 
integral approximation, to guarantee approximation over &/, or equivalently 
over J2x(I 0)OO ) ■ 

4. Discussion on approximation 

We particularize in this section a proposal for approximation, and we 
analyze the main properties of the error in the time and frequency domains. 
These properties will be still valid for any other approximation in the graph 
topology over 

4-1. Proposal for a constructive approximation 

We propose in this subsection an effective and constructive approximation 
using elements in Jtf s (l$ 1 ^ 2 ). We start this proposal by the particular case 
of 0„(O- We denote C k n = F ^ yT , MO = ^(-a-UnO, for £ e]0,l], and 
^o(O) = 0. 

Lemma 4. Consider the sequence #o,n(") in J^(Io^) described by 

e 0l n(t) -E C ^° (") e- ak \l-e- at ) n - k , t G V 
fc=o 

and elsewhere. Then 0o,n( - ) uniformly converges to 9q(-) for the stf -norm. 
Proof. Define 

MQ) = o (-a- 1 \ng), g e]0,l], (28) 

and V'o(O) = 0. The function Vo(") nas a bounded step discontinuity for 
g = e~ ai? . Since e _m? is irrational, there exists k in 0, 1, . . . , n — 1 such that 
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n/n < e Q1? < (k + l)/n. We take ipo,c(') an y continuous function over [0, 1] 
satisfying 

tpoM = ^o(^), Q e [0, «/n] U [(k + l)/n, 1], 

that will remove this discontinuity and will be close as desired to ipo{') f° r 
the ^fi(Io,i)-norm. In particular ipo,c(') w iH satisfy 

ipo,c( k / n ) = i>o{k/n), k = G,l,...,n. 

We now approximate as close as we please the function ipo,c (or equivalently 
i/jq) by Bernstein polynomials 

^0,app(f?) =J2 C n^ (-) Q k (l - QY~ k - 
k=0 ^ ' 

By the inverse transform of ( |28j) . we obtain 

VW=^o,a PP (e- at ), 

which is in J^(Io,#) since ^o(O) — 0. Furthermore, from Lemma El we know 
that #o,n(") converges uniformly, by construction, to #o( - )- '-' 

Previous lemma gives an approximation of 8q in 3&(9o,e n ), where positive 
upper bound e n can be reduced arbitrarily when n increases. This algorithm 
can be easily generalized to approximations in J^(I 0jOO ) or of the form (127|) . 
by modifying the domain of definition for (12 81) . For the more general case of 
&\{-), we propose the following constructive solution. 

Lemma 5. The sequence of functions defined over J%,(Io t &) by 

<M*) = £ c » *a (£) (e- - e— )* (1 - B-)-*, 

^ ' fc=0 v ' 

for t G Io,tf> anc ^ 

$ A (/i) = ^(-a- 1 ln((l - e-«*)/x + e"^)), // G [0, 1], 
uniformly converges to #>(•). 



19 



Proof. From Lemma [2J define 



$ A (» = 9 x {-a' 1 ln((l - e-^/i + e""*)) 

with /i G [0,1]. Polynomial approximation $A,app(') can be obtained from 
Bernstein polynomials, that is 

$A,app(//) = C n**( ~ ) A 1 " ^ n " fc - ( 29 ) 

The inverse transformation for $A,a PP yields to 

M*> = 7^— h-r E c «*^ (s) (e "°' " e_<rf) ' f 1 - 6 ""')""'- 

which uniformly converges to 9\(-). □ 

This approximation writes like a sum of elements in J^(Iq^) and 6*o(')- With 
Lemma H] and approximation of (-), such an approximation is defined over 
J^(Io,#), and writes like in ( 12"5|) . Some basic considerations can be made 
on the order of this approximation. We know from [21] or [22] that, given 
e > 0, there exists r\ such that for any G [0, 1], — yu 2 | < rj implies 

|$aOi) - < $g, and that 

|$A(/i)-<fA, a pp(^)|<^+ ll<1,Xl|/ ^ 



2$ 2$r] 2 n 

There exists /3 a positive bounded real number, such that 



<P\m- 



2M 



/i 2 |, with for instance f3 > 1 e _ e jJ^ . This in turn implies that rj < ^e 
Taking the maximal admissible value for rj, we finally obtain that the ap- 
proximation order n satisfies 

n > 

e 6 

From this simple consideration coming from the use of Bernstein polynomials, 
we obtain, to guarantee a norm upper bound e for the error, a condition on 
the order of the approximation. This condition is however quite conservative, 
and in practice, the order can be chosen iteratively, as illustrated in Fig. [3] 
and Fig. gj 
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Figure 3: Kernel approximations #i ia pp(i) using exponentials in the time domain of the 
kernel 9\{t) (continuous line), for orders n — 5 (dashed) and n = 10 (dot-dashed). 
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0.2 - 



1 23456789 10 
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Figure 4: Representation of the norm error \\9\ — 9i^ n \\s^ with respect to the order n of 
the approximation, for a — 1 (star) and a = 1/5 (circle). 

Properties in time and frequency domains 

Let / app be a kernel approximation in £%(f,e), for a given e > and 
/ e Jf(I# u # 2 ). The corresponding output is y &pp {t) = (/ app From flBJ), 
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for any 1 < p < oo and u G J^,(R + ), 

\\y ~ l/applU, < (3°) 

In other words, the output error e(t) = y(t) — w app (t) can be made arbitrarily 
small for the Jz^,-norm by a suitable choice of the arbitrary bound e. Such a 
property includes the case of persistent inputs in «if 1 (R + ). For the particular 
case p = 2, we verify that the proposed approximation holds also in the 
frequency domain, since 

sup |/(5) - / app (s)| = 11/ - f app \U x < e. (31) 

Res>0 

Theorem 4. Let f £ be a given distribution, and / app an approximation 
in 33(f,e), for some small e > 0. Then 

11/ - /applUoc = SU P \f(j U ) ~ /appCMI < £, 

and /or a// u in R, 

|arg(/0' w )) ~ ar g(/ap P (jw))l < e- 
Proof. Since / app e M(f,e), and 

11/ — fappWjifoc — ||/ — /app Hi/; 

module inequality is trivial. To show that approximation holds also for the 
phase angle, we denote / app = / + e app . Hence 



\f(ju)\ |1 - e^)-^(-»| - \e app (ju)\ 



<e, 



where (p(u) = arg(/(ja;)) and y? ap p( w ) = arg(/ app (ju;)). This implies that 

|/0'w)| |l-e* (v(w) ~ ,p "« ,(w)) | 

can be made as small as we please for all u, which in turn implies that (p(u) 
and (p a pp(uj) are arbitrarily close for the Jzfoo-norm. □ 

This frequency property is illustrated in Fig. |5l where an approximation with 
order n = 5 in J^(I 0j i) of 6\{joS) is plotted. The approximation of distributed 
delay in the graph topology over srf yields to an approximation in both time 
and frequency domains. This is a strong property, that turns to be central 
in control problems. 



22 



20 




Frequency w 




10° 10 1 10 2 10 3 



Frequency w 



Figure 5: Bode diagram of Oi(jto) and its kernel approximation 9i tapp (ju}), with order 
n = 5. 



5. Application to control of time-delay systems 

Distributed delays appear naturally in the control of time-delay systems. 
We consider here as control application, the stabilization problem, and we 
illustrate the approximation method to realize the distributed time operator. 

5.1. Stabilization 

In the stabilization problem, we determine a control law of a given plant, 
such that the closed-loop system is stable. For this, consider a distribution p 
with a coprime factorization (n, d) over and define the closed-loop system 
described in Fig. [6j with a compensator c defined over the quotient field of 





c 






V 









2/2 



Figure 6: Feedback System. 
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We assume here that the plant p admits such a coprime factorization over srf ' . 
This is an assumption, but most dynamical systems fulfill it. The dynamical 
equations of this system are 



H(p,c) 



"2 



H(p,c) 



l+pc 

PC 



—pc ■ 
l+pc 
V 



. l+pc l+pc. 



(32) 



We say that c stabilizes p, or the pair (p, c) is stable, if the matrix H(p, c) G 
£/ 2x2 . All internal signals in the closed-loop are bounded for any bounded 
exogeneous inputs U\ and w 2 . A necessary and sufficient condition for (p, c) 
to be stable is the existence of some n c and d c in A such that 



n * n c + d * d c 



v. 



(33) 



with v a unit in stf . Since (n, d) is a coprime factorization, (|33|) holds, so 
we conclude that a stabilizing compensator c has a coprime factorization 
(n c , d c ). A coprime factorization (n c , d c ) over ^ of a stabilizing compensator 
for p, includes, in general, some distributed delays. Approximating such a 
compensator as described in Section [3] yields to the approximated controller 



c app with factorization (n Cjapp , 4 >app ), where n Ciapp G 



and rf Cjapp G 



38(d c ,Ed), for some given e n and Applying this approximation in the 
control, we would know the conditions on c app , such that the pair (p, c app ) 
is stable. The positive answer comes from a direct application of small gain 
theorem for BIBO-stability. 

Lemma 6. The pair (p, c app ) is stable if 

m&x(e n ,e d ) < \\(n d)\\~J . 

Proof. Let (n c , d c ) be a coprime factorization over srf of a stabilizing con- 
troller. From f )33|) . we have 

n * n c + d * d c = 5. 



The approximated controller c app with factorization (n c app , d CAPP ) yields in 
closed-loop 



n * n Cjapp + d*d, 



c,app 



[n d) 



^c ^-c,app 
dc ^c,app 



(34) 
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The approximated controller will stabilize the plant p if and only if the right 
hand side in (j34p is a unit over srf . Since s$ is a Banach algebra, a sufficient 
condition is 



(n d) 



dc d ca pp 



< 1. 



Via approximations, n capp e 3&(n c , e n ) and c/ capp G 3§{d c ,e^), so that the 
above inequality yields the sufficient condition 

|| (n d)||^max(e n ,£ d ) < 1. 

□ 

As a first comment, remark that it is always possible to determine an ap- 
proximation c app such that (p, c app ) is stable. The counterpart will be in 
the order of approximation, that will increase when the required accuracy 
vanishes. Small gain theorem, which is still valid over any Banach algebra, 
gives us a sufficient condition on the approximation accuracy to guarantee 
robust stability. This condition helps us to determine the order of the ap- 



proximation. Note that from [31[ where the conservativeness of small gain 
theorem for BIBO-stability was studied, a converse statement for Lemma [6] 
holds, like in the case. This highlights the weak conservation of such 
a condition for robust stabilization. As an application, consider the plant 
y(s) = p(s)u(s) given by 

m = (35) 

A coprime factorization writes h = jtj, d = ^j-, since 

n(s) 2 e 1 + d{s) (l + 2 O^s)) = 1. 

Hence, a stabilizing compensator for ( |35l) is 

u(t) = -2(6 1 *u)(t) + 2e 1 y(t). 

From Lemma El a sufficient condition for robust stability with approximation 
of 6*i is that Ed < ^rg- From Fig. HJ we see that a first order may be sufficient 
for stability purpose. In practice, a sufficient accuracy is obtained for a 
5th order approximation over J^(I 0j i), as can be seen in Fig. [31 where step 
responses are plotted. Note also that if we modify the coprime factorization 
(n,d), we could obtain lower order approximations. 
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Figure 7: Step responses for the desired closed-loop (continuous line) and for the closed- 
loop system with approximation of distributed delay (5th order, dash). 

6. Conclusion 

In this paper, we proposed a general framework for rational approxima- 
tion of distributed delay operators. The properties of this approximation 
were studied in both time and frequency domains. The effectiveness of this 
method was shown in simulation on the stabilization and the finite spectrum 
assignment problems, for general linear time-delay systems. We expect that 
this work will bring some new light in the understanding of distributed delay 
approximation, and more generally in approximation for the control of linear 
infinite dimensional systems. Such an approximation will provide founda- 
tions for a renewal of interest in control synthesis with distributed delays. 
Distributed delays inherit from integral control good robustness properties, 
and are at the core of numerous results in robustness analysis, optimization 
and control for time-delay systems. 
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